arXiv:1505.05237v2 [hep-ph] 27 Oct 2015 


High-energy electrons from the muon decay in orbit: radiative 

corrections 

Robert Szafron, Andrzej Czarnecki 
Department of Physics, University of Alberta, Edmonton, 

Alberta, Canada T6G 2G7 

October 28, 2015 


Abstract 

We determine the 0{a) correction to the energy spectrum of electrons produced in the decay of 
muons bound in atoms. We focus on the high-energy end of the spectrum that constitutes a back¬ 
ground for the muon-electron conversion and will be precisely measured by the upcoming experiments 
Mu2e and COMET. The correction suppresses the background by about 15%. 


In matter, muons decay differently from antimuons. Although the decay rates are very similar [1], 
negatively charged y~ can bind with nuclei. The nucleus exchanges photons with the muon and the 
daughter electron, rearranging the energy distribution. In this paper we find how this rearrangement 
is affected by the real radiation and self-interaction on the muon-electron line. We predict the energy 
spectrum of the highest-energy electrons, interesting both theoretically and experimentally. 

For a theorist, the muon decay is the simplest example with which to understand the gamut of binding 
effects, including the motion in the initial state, interplay of the binding and the self-interaction, and 
the recoil of the nucleus. Experimenters have recently studied the bound muon decay (decay in orbit, 
DIO) [2] with a precision sufficient to probe radiative corrections, later evaluated in [3]; however, these 
studies concern only the lower half of the spectrum, largely accessible also to a free muon. 

Interestingly, the energy range of electrons produced in the DIO reaches to about twice the maximum 
possible in a free muon decay. When the muon decays in vacuum, momentum conservation requires that 
at least half of the energy be carried away by the neutrinos. In the DIO, the nucleus can absorb the 
momentum without taking much energy because it is so heavy. 

The high-energy part is important for the upcoming searches for the ultra-rare neutrinoless muon- 
electron conversion, COMET in J-PARC [4] and Mu2e in Fermilab [5]. Designed for a sensitivity better 
than one exotic conversion in 10 16 ordinary muon decays, they will collect large samples of events with 
high-energy electrons. A reliably predicted spectrum is needed to distinguish the exotic signal - an excess 
of electrons at maximum energy - from the Standard Model background. 

Predicting the DIO spectrum is a challenge because both the decaying muon and the daughter 
electron interact with the Coulomb field of the nucleus. A numerical calculation with Coulomb-Dirac 
wave functions is possible [6] provided that self-interactions (photons attached to the muon and the 
electron) are neglected. How can they be included? In the lower half of the spectrum the muon and 
the electron can be treated as nearly free and the binding effects can be factorized. Then the radiative 
corrections, known for a free muon, are convoluted with a shape function that parametrizes the Coulomb 
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field effect [3, 7]. Here we construct an expansion around the end-point and employ it to find radiative 
corrections also to the high-energy part of the spectrum. 

Accounting for the external Coulomb field in charged-particle propagators is called the Furry picture 
[ 8 ]. In this formulation, and still ignoring radiative corrections, a single diagram, shown in Fig. 1, 
describes the DIO. We shall demonstrate that the bound-state radiative corrections are easiest to evaluate 




Nucleus Furry picture 

Figure 1: Muon decay in orbit (DIO). Dashed lines denote Coulomb photons exchanged between charged 
leptons and the nucleus. The right panel shows the same physics using double lines for charged leptons 
propagating in the Coulomb field. 


near the high-energy end of the spectrum, the most important part for the new experiments. For now we 
neglect the nuclear recoil and structure, and treat the nucleus as an inhnitely-heavy point source of the 
Coulomb field. We denote the electron energy with E; its maximum value is E max ~ ^1 — ^, 

where is muon mass, Z is the atomic number, and a ~ 1/137 is the fine-structure constant. The 
DIO spectrum near its end-point can be expanded in the small parameter A = Em * x-E , 


dr 
r 0 d E 


Y,Bi ik K{*ZaY Q k , 


( 1 ) 


G 2 m 5 

where r 0 = is the free-muon decay rate and G P is the Fermi constant [9, 10]. Powers of 7 xZa 

parameterize photon exchanges with the nucleus and a/- 7 r arises from radiative corrections on the charged- 
lepton line and the vacuum polarization. The first non-vanishing term has i — j = 5 and k = 0, with 
H 550 = ~ 0.21. Higher order coefficients B may have logarithms of Za and A. 

Corrections to this leading behavior have several sources. The large momentum transfer to the 
nucleus probes its interior. The finite nuclear size, already included in [ 6 ], causes the largest correction. 
We will comment at the end of this paper on how to include it in our formalism. The finite nuclear mass 
introduces a recoil effect, also evaluated in [ 6 ]. It affects the coefficients B only slightly but it shifts the 
end-point energy E max . 

We shall exploit a theoretical similarity between the DIO and the photoelectric effect to control 
higher-order binding effects. They generate powers of ixZa [11, 12] rather than Za. Indeed, a nu¬ 
merical evaluation for a point nucleus with Z = 13 (as in aluminum, the planned target in COMET 
and Mu2e) finds a —21% correction, consistent with 137ra = 0.3. Logarithmic enhancement starts with 
(ttZ a) 7 ln(Za). Fortunately, these large effects, slightly suppressed by the finite nucleus size, are summed 
up in the numerical evaluation [ 6 ]. 

Finally, the most challenging corrections result from radiative effects that are the subject of this study. 
Before delving into the physics of the end-point, we present our main result. Close to the end-point, 
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including radiative corrections, the DIO spectrum for aluminum is 


— « 1.24(3) x 1(T 4 x A 5 ' 023 . (2) 

r 0 dE w w 


To illustrate the importance of the new corrections we consider the last 150 keV of the spectrum (the 
typical planned resolution of Mu2e and COMET). Radiative corrections reduce the number of events in 
this bin by 15%, a welcome reduction of the background, comparable in size with higher-order binding 
effects. 

In the remainder we explain the origin of such a large effect. We begin with the tree-level behaviour, 
appropriately expanding the lepton wave functions. We find that an exchange of a single, highly virtual 
photon gives the electron an energy of the full muon mass. 

The relativistic electron is described by a plane wave distorted by the Coulomb potential V; to the 
first order, 


% (q) = u(p) 


S 3 {p- q) + f ({p~ qf) 


( 3 ) 


where u(p) is a spinor solution of a free Dirac equation and the four-potential in momentum space reads 


V 



( 4 ) 


A muon bound to a nucleus with Z <C 137 is nonrelativistic. Nevertheless, we will need the first 
relativistic correction to its wave function, just like in the classic analysis of the photoelectric effect [13], 


= i’NR(q) ^1 + U ( P )> ( 5 ) 

where ip NK (q) = nonre l a tivistic momentum-space wave function of the IS ground 

( Zam \ 3 / 2 

state with \H(0) = ( ) ; u(P) is the four-spinor of a muon at rest, P = (m M , 0). 

We now consider separately the contributions of the two terms in the electron wave function (3). The 
delta function term forces the muon momentum in (5) to be large, q = p ~ Thus we neglect Zam M 
in the denominator of i/> NR and find 


4>(q) 


(2tt) 3 T(0)- 


'+ 4-m. 


-f (q 2 ) u(P). 


( 6 ) 


This is visualized in Fig. 2a: the muon, before decaying, transfers momentum q ~ to the nucleus 
through a hard space-like photon. It is here that the relativistic correction to the muon wave function is 
important. 

The second term in (3) refers to an electron scattered on the nucleus. Now the muon momentum, 
not restricted to large values, has its typical bound-state size q ~ Zam^, negligible in comparison with 
p ~ m M . We use lim^o = (2n) 3 5 3 ( q) to approximate the muon wave function, 

^(q)^(27r) 3 ^(0)8 3 (q)u(P). (7) 


This is shown in Fig. 2b, where the hard photon is exchanged after the decay. 

The two diagrams in Fig. 2 add up to the leading contribution -B 550 in (1). In both cases any energy 
unused by the electron (~ A) is taken up by the neutrinos and not transferred to the nucleus. Counting 
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Figure 2: Furry diagram expanded in Za. Crossed circles indicate insertions of the weak interaction 
transforming the muon into an electron; the emitted neutrinos are not shown. These two amplitudes 
give rise to the highest-energy electrons. 


neutrino momenta in the integrated matrix element explains the leading energy dependence in (1), 

r d_u d_^o_ s ( m A _ v 0 -V 0 ).. .f.. ,~j ~ A 5 . (8) 

J v0 Vq 

Having understood that only two diagrams describe the end-point behavior, we are now ready to 
evaluate radiative corrections. In the Furry picture there are two groups of virtual corrections, shown in 
Fig. 3. We expand them in Za just like the tree-level diagrams, but in addition to wave functions (3, 5), 
we need also the Coulomb-Dirac Green’s function [14], 


( E-,p,p ') 


S 3 

f — m 

+j— V (( p-p 0 2 ) -j— 

® — m V / © — 


The expansion (9) reduces radiative corrections in an external field to a set of loop diagrams that we 



evaluate analytically [15]. This approach can be extended to higher-order corrections. 

The diagram in Fig. 3a adds the Uehling correction to the Coulomb potential [16] and modifies the 
photon propagators in Fig. 2. In muonic aluminum the range of the Uehling potential exceeds the Bohr 
radius. It strengthens the attractive force and increases the muon wave function at the origin, 

’F(O) -»• ’F(O) (l + ^<5 0 ) . (10) 

For aluminum we find 5 o = 3.27. This correction reflects the running of the coupling a up to the average 
muon momentum scale m^Za. 

Vacuum polarization loops on the highly virtual photon propagators are related to the running of a 
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up to the hard scale m M . They enhance the tree-level decay rate by a factor 1 + -<5 V p, with 

5vp = +0.12 « 6.1, (11) 

3 m e 9 

where the term 0.12 arises from a muon loop. 

Another correction comes from the real radiation. Diagrams represented by Fig. 4 are expanded in 



Figure 4: Furry diagram for the real radiation correction. 


the same way as virtual corrections, using (9). Near the end-point the eikonal approximation suffices; 
by energy conservation the real photons must be soft, 0 < < m M A. 

The sum of virtual and real radiation is finite, 

■ 75 — = <5 H + i$s In A, (12) 

-0.550 

where the vacuum polarization corrections (5 V p and 5 0 are included with hard self-interaction effects in 
<f H = 6.31 — Ulnyfru and S s = 21 n^hi — 2 is a soft correction. The latter can be exponentiated [17] 
(similarly to the free-muon decay [18]) and vanishes when A —> 0, 

B 550 H--B 551 —► -B 550 (a^ 5 ® H— Sh) , (13) 

instead of unphysically diverging as In A. It increases the exponent of A and suppresses DIO events near 
the end-point. The relative decrease is inversely correlated with the energy resolution: the number of 
electrons in the end-point bin of 1 (0.1) MeV is reduced by 11% (16%). 

The final-state electron is relativistic, E m e , so its structure function [19, 20] is insensitive to 
Coulomb corrections. A convolution with the free-decay spectrum confirms the coefficient — of In 
due to collinear photons. Together with the vacuum polarization in (11), this explains the logarithmic 
part of the hard correction. 

That log is largely cancelled in the sum with the wave function correction in (10) and <5 H = —2.9 
reduces the end-region by only a fraction of a per cent. We thus neglect the unknown hard corrections 
O ((a/-7r) 2 ) in the error estimate, dominated by the nuclear-size effects, discussed below. 

There are now two complementary studies of the end-point spectrum. Here, we have computed 
radiative corrections (RC) assuming a point nucleus and considering only the one-Coulomb exchange. 
Ref. [ 6 ] did not have the RC but included the nucleus structure, recoil, and multiple Coulomb interactions. 

In order to combine these results, we observe that the most important - soft - correction is universal, 
not sensitive to any interactions with the nucleus. The hard correction is tiny, so treating it also as 
universal is well within our final error estimate. 

In the discussion of the uncertainty we specialize to aluminum but the discussion can be applied to 
other nuclei, so we keep the Z dependence explicit. We assume a Fermi charge distribution, 


Q = 


Po 

1 + exp ^ 
x ao 


(14) 
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with ao = 0.569 fm and ?’o = 2.84(5) fm [21]. 

The finite size affects the nucleus form-factor, defined as a ratio of Fourier transforms of potentials 
from the extended (14) and the point-like (4) charge distributions, 

V P (k 2 ) 

F p (k 2 ) = — 7 —r-- > 0.64 for k 2 = m 2 . (15) 

V (k 2 J 


The DIO spectrum for a finite nucleus has an expansion analogous to (1), but with coefficients that 
depend on the density p. Its leading term near the end-point [ 6 ] is 


OO g 

B£ j0 \F p (ml) nZa\ j =8.98 xl0~ 17 (^) . (16) 

1=5 


This result includes exchanges of many Coulomb photons, in addition to the single hard exchange to 
which we have found the radiative correction. We estimate the magnitude of the multi-Coulomb part as 
a fraction / = F p (m^) ttZ a ~ 0.2 of (16). 

Hard radiative corrections to this part are missing. To be conservative, we are not assuming that 
they involve a cancellation that has suppressed <5 H - Corrections on the order of the collinear logarithm 
translate into a relative error of about ^|^/ln ^ ~ 0.7%. In addition, experimental errors in the charge 
distribution parameters (14) introduce a 2% uncertainty [ 6 ]. Summing them in quadrature, together with 
the sensitivity to the scale involved in the exponentiation of soft effects, we arrive at an error around 
2.5% in the end-point spectrum. 

The result (16), multiplied by the new correction (13), leads to our prediction for the end-point 
spectrum, ( 2 ). 

To summarize, we have determined the correction to the high-energy tail of the DIO energy distri¬ 
bution and its remaining uncertainty. Key to this improvement has been the simplicity of the leading 
amplitudes that turn out to arise from a small number of hard-photon exchanges. This line of reason¬ 
ing can be extended to higher-order binding effects, at least for a point nucleus. For a realistic charge 
distribution, a numerical evaluation of loop diagrams will be necessary. However, the leading radiative 
correction has now been established with good precision. Its sizeable negative effect on the DIO will 
make any observed event near the end-point an even more convincing signal of New Physics, a discovery 
we eagerly anticipate. 
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